function [ H1 ] = ApplyHConservaty( H1, R )
    cH = sum(H1);
    cH = repmat(cH, R, 1);
    H1 = H1 ./ cH;
end

